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mean curvature in Riemannian submersions fibered by flow lines of a 
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1 Introduction 

Recent papers devoted to the study of CMC surfaces in certain homogeneous 
three-manifolds are based in the description of these ambient spaces as Rie- 
mannian submersions over constant curvature model surfaces. For instance, 
this is the case of pQ, [3], [7] and [8]. In particular, in [3] the authors obtained 
CMC graphs in the Heisenberg space regarding it as a submersion over M 2 
fibered by geodesic flow lines of a Killing vector field. The goal in these works 
is to extend classical results about CMC surfaces in Euclidean space as well 
as more recent results in nonflat space forms to a more general setting. 

One of the main issues in developing a theory for CMC hypersurfaces in 
general Riemannian ambients is the existence of examples. Methods which 
rely mainly on geometric constructions may fail if the ambient space lacks 
appropriate symmetries or structures. However, the problem may be solvable 
once it is reformulated in analytical terms as the existence of CMC graphs 
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for a suitable notion of graph. This is the case of Riemannian manifolds 
carrying a Killing vector field where the natural notion of Killing graph has 
been defined under additional assumptions. 

The Dirichlet problem for prescribed mean curvature Killing graphs in 
ambient spaces endowed with a Killing field with integrable orthogonal dis- 
tribution was first solved for CMC surfaces in [5]. Then, it was extended 
in [5] to hypersurfaces with prescribed mean curvature function. Under the 
integrability assumption, the ambient manifold has a warped product struc- 
ture with one of the factors giving rise to a totally geodesic hypersurface 
foliation. 

In this paper, we consider a generalization of [3] and [5] to Riemannian 
submersion ir : M n+1 — > M n whose vertical fibers are given by flow lines of a 
Killing field. Thus, the normal distribution to the Killing field may fail to be 
integrable. Our aim is to show that a natural setting of the Dirichlet problem 
for Killing graphs (defined in Section 2) with prescribed mean curvature 
function in this context is to consider these as leaves transversal to a solid 
cylinder of the flow lines that project on a compact domain on the base 
of the submersion. Using this approach, we give a unified proof of known 
and completely new existence results in a wide range of ambient Riemannian 
manifolds. Among the ambients for which this paper applies, we should 
mention higher-dimensional Heisenberg spaces and odd-dimensional spheres 
submersed in complex projective spaces. 

The existence part of our result is proved using the continuity method 
for quasilinear elliptic PDE. In order to obtain apriori estimates essential to 
this method we use Killing cylinders as barriers. Given a domain Q in M 
with compact closure and boundary T, the Killing cylinders over T and 
are, respectively, the subsets K = 7r _1 (T) and M = ir^^Cl). We denote by 
if cy i the inward mean curvature of K and by Ric^ the Ricci tensor of M. 

With the above notations we have the following result. 

Theorem 1 . Let Q C M be a domain with compact closure and C 2,a bound- 
ary. Suppose that H cy \ > and mi^RicM > — n infr if,? yl . Let H G C a (fl) 
and G C 2,a (T) be given functions and l: Q — ► M C M be a C 2,a immer- 
sion transversal to the vertical fibers such that it o i = id\^. If 

sup | if | < inf if C yi, 
n r 

then there exists a unique function u G C 2,a (Cl) satisfying u\y = <p whose 
Killing graph S has mean curvature H. 
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The hypothesis on the existence of an immersion i is used simultaneously 
to introduce a set of coordinates well suited to the problem and to define 
properly the notion of Killing graph. In terms of these coordinates, it may 
be rendered evident that the ambient metric is stationary. Moreover, is 
used as barrier to producing an initial minimal graph by the direct method in 
Calculus of Variations. In higher- dimensional Heisenberg spaces there exists 
a minimal leaf transverse to the flow lines of the vertical vector field. Thus, 
in this particular case there is no need of the hypothesis. By contrast, if we 
consider the example of odd-dimensional spheres submersed in the complex 
projective spaces, it is not guaranteed that always exist such minimal graphs 
with respect to the Hopf fibers. 

We remark that submersions with totally geodesic fibers constitute an 
important example where we may construct initial Killing graphs. In fact, if 
we also assume that the Killing cylinder Mq over Q is geodesically complete, 
then geodesic cones with boundary in K and vertex at the mean convex side 
of K may be taken, after smoothing around the vertex, as initial Killing 
graphs. Thus, we may rule out the hypothesis in this case. 

This paper is organized as follows. In Section 2, we fix notation and made 
precise the notion of Killing graph. We deduce the mean curvature equation 
and define adapted and basic reference frames crucial in the subsequent anal- 
ysis. In Section 3, we present some basic geometry of Killing cylinders. In 
Sections 4 and 5 we construct analytical barriers to obtain height and bound- 
ary gradient estimates. Section 6 is devoted to the proof of interior gradient 
estimates based in the technique of normal perturbation of the graph due 
to Korevaar [13]. The continuity method and the existence of the minimal 
initial solution are presented in the final section. 

2 Killing graphs 

Let 7r: M n+l — > M n be a Riemannian submersion such that the leaves of 
the vertical foliation are the trajectories of a nonsingular Killing vector field 
denoted by Y e £(M) . Let Q C M be a C 2,a domain with compact closure. 
We assume that the integral curves of Y in 

M := 7r -1 (fi) 

are complete lines. Since the hypersurfaces we work with are graphs over O 
along the integral curves, when these curves are circles we may pass to the 
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universal cover of M without loss of generality. 

Let i: Cl — > M be an immersion satisfying not — id^ such that the 
hypersurface So = l(JT) is transversal to the flow lines. The initial values for 
the flow ^: H. x £ — > M of Y are taken at E , i.e., S corresponds to the 
level hypersurface s = for the flow parameter s. Set ^ s = \&(s, • ). Then, 
the level hypersurfaces S s = X I / S (S ) constitute a foliation of M by isometric 
hypersurfaces. 

Fix a local reference frame vi, . . . , v n on Cl and set 

(Tij = (Vi.Vj). 

Let . . . , v n be the corresponding local frame on S , i.e., Vj(p) = t*Vj(x) 
if x G f2 and p = By means of the flux \I/ we define a local frame at 

q = * a (p) in M by 

9 a (g) = ^ tt(s,p) = r(*(s,p)) = #*(s,p)d s (p) 

and 

= (* a o t)vi(x). 

Let Di, . . . , D n in M denote the basic vector fields 7r-related to vi, . . . , v n . 
If q — ^(s,p) for p G S , then ir(q) = n o ^(s,p) = n(p). Therefore, 

A(g) = *.(s,p)A(p) 

since is horizontal and 

7r*(g)#*(s,p)A(p) = (tt o tf)»(s,p)A(p) = 7T*(p)A(p)- 
That 7r is a Riemannian submersion yields 

Setting 

A,:=/ 1/2 9 s , 

we complete a local reference frame Dq, D\, . . . , D n on M where / := 1/|^| 2 
does not depend on s since Y is a Killing field.. 

We extend the frame v 1? . . . , v„ adapted to the leaves S s to a frame 
Vs, . . . , v„ in M by adding the gradient vector field Vs of the function s. 
Using 

7r*(?)vi = 7r*(p)i*Vj(a;) = Vi(x) = -K*{q)Di 
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and 

l = d s s = (Vs,d s ) = f- 1 / 2 (Vs,D ), 
we have that the two frames considered on M are related by 

Vs = f 1 l 2 D + a^D j (s)D l 
\7j = 6iD + A- 

The functions 5« are independent of s since 

5, = (v<(g), A,(g)> = A>(p)> = (Vi(p), A)(p)>- 

Thus, from 

= v,( S ) = (V S ,v J )=/ 1 / 2 5,+ J D J ( S ) 
we conclude that the functions Dj(s) are also independent of s. 

The Killing graph £ = £ u of a function u G C 2 (0) is the hypersurface 

S u = {fW P ),p):pGS }, 

where u is seen as a function on S by taking u(p) = u(x) when 7r(p) = x. 
Since £ can also be considered as given by the immersion 

l u : x G f2 I— > ^(^(x), l(x)), 

its tangent bundle is spanned by the vector fields 

(i u )*Vi = Vi(u) ^ s + (* o t ) +Vi = Vi(u) <9 S + Vi = r 1/2 Vi{u)D + Vi. (1) 

We may regard function in M by means of the extension 

u(q) = u(x) if ir(q) = x. (2) 

Thus Dq{u) = f l l 2 d s u = 0, and hence 

AO) = Vi(u) - SiD (u) = Vj(w) = Vj(tt). 

Therefore, we have using ([T]) that 

(0*v l = (r 1/2 A(«) + 5 l )^o + A. 
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It follows easily that a unit normal vector field to £ pointing upwards is 

N = ^(f 1/2 D -uW j ), (3) 

where 

u j := a ij Di(u - s) (4) 

and 

W 2 :=f + a ij u i u i = f + u% 

for Ui := GijV? . Notice that u? and W can also be seen as functions on M 
since they are independent of s. 

2.1 The mean curvature equation 

To compute the mean curvature of £ assume for simplicity that the tan- 
gent frame Vx,...,v n is orthonormal at x G Q. Hence, the basic frame 
D , Dx, . . . , D n is orthonormal at points of 7r~ 1 (x). Thus, 

fl/2 iff I 

(V Do N,D ) = ±— {VdoDqjDo) _ w{Vdo D v D ) = —(V^Da^Dj), 

where V denotes the Riemannian connection on M. We consider on M the 
vector field 

Du := u 3 Mj = a l3 Di(u — s)vj. 
Since Vd Dq is a horizontal vector field and 7r*(u J Dj) = u J Vj, we obtain 

(V Do N,D ) = ^(<k*V Do D ,Du). 

By the well-known O'Neill submersion formula [9j, we have 

Vn.D^iV^D^ + ^Dtf. (5) 

Thus, we obtain for k > 1 that 

f 1/2 - - /u J \ Du 

(V Dk N,D k ) = - J w (V Dk D k ,D ) - (V Dk [ w D 3 ),D k ) = -(V Vfc — ,v fc ) 

where V denotes the Riemannian connection on M. We conclude that 

Du 1 

nH = -divjv/ N = div M -jrr - — (^V^ A), Du). (6) 
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Denote the covariant derivative in M of Du = vpMj by 

and set ■= ajkU 1 .^ Computing at any point the divergence in (jHJ) gives 

i- Du ik/T-, Du \ 
div M — = a (V v — ,v*> 



a ik 

W 2 



a ik ■ 1 

— {W 2 {u ] . l y j ,M k ) - -v<(/ + u l ui)u>a jh ) 



On the other hand, 



Vi(/) = -/ 2 v,(r,r) = -/ 2 (5 lJ D (r,r) + A(^^)) 

= -2/ 2 (V Ci Y,Y) = 2/ 2 (VyY, A) = 2f{V Do D , A) 
= 2/(7r»V A) £>o,v i ). 

Thus, the mean curvature equation ([6]) becomes 

- (/ + ^ 2 ) (7r*Vc A, Du) = B, (7) 

where 

A ik := W V fc - and B := n# W 3 . 

We define the operator 

Q[u] = - (/ + ^ 2 )(vr,V A) A,^)). 

Therefore, we have shown that £ is a hypersurface with prescribed mean 
curvature function H(x) and boundary condition if u is a solution to 
the Dirichlet problem 

Q\u] = nH 
I ^ (8) 

where V = dQ. The boundary of £ is the Killing graph over T of 0. 
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2.2 A commutation formula 

In this subsection, we give a commutation formula for second covariant 
derivatives that allow to conclude the ellipticity of the quasilinear operator 
defined in the proceeding one. 

Since Du = 7r*V(tt — s) from (j3J we obtain using (jSJ) that 

V Vk Du = V Vfc 7r*V(-?i - s) 

= Tr*(V Dk y(u-s)- V Dk (V(u-s) 1 D )D ) 

= a il 7r*((V Dk V(u - s), A) A - (V(u - s), A) (V Dfc A, A) A) 

= ^(V^, A (u - s) + (£>„(«) - D (s))(V Dk D h A)) ?r*A 



^(V^a («- a) - ^([A, A], A» v,. 



Therefore, 



= = ^(V^ D; (u- s) - ~/ 1/2 <[A, A], A)) 

= V^ >0 . (« - s) + - ljk (9) 

where j k j := /^([Aj A?']> A)) is skew-symmetric. Since the Hessian is 
symmetric we conclude that 

Uj-k - u k .j = j jk . (10) 
Under the convention for u established in ([2]) we use the standard notation 
Ui = Di(u), u % = u^Uj and u^j — (V^Vw, Dj). 
Then the matrix (ui-j) is related with the Hessian matrices (ui-j) and (sj ;J ) by 

Ui;j = (V Dj Vm, A) - (V Dj Vs, A) + ^7ij = ~ + ^ Tir 

Hence, the principal part of the mean curvature equation (jSJ) is given by the 
matrix {A 1 ^). This matrix is positive-definite. Indeed, we have that 

f\i\ 2 < A* fa < W 2 \i\ 2 . (12) 
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3 Killing cylinders 



The Killing cylinder over T is the hypersurface K = n 1 (r). Thus 

K = {V(s,i(x)) : s e R, x e r} 

is ruled by the flow lines of Y through t(r) C S . 

We denote by fj the inward pointing unit vector field normal to K. Clearly, 
fj is a basic vector field and 7r* fj = 77 is the unit normal vector field to T in M 
pointing inward. We work with a tangent frame satisfying that vi = rj and 
v 2 , . . . , \i n are orthogonal to v x . In particular, their horizontal lifts Di verify 
along K that D x = fj and Dj, 2 < j < n, is tangent to K. Set 

f(V Y Y,f}) = (V Do D ,fj) = K, 

where k can be seen as a function on Cl. In fact, 

f(V Y Y,r ] } = -f(V n Y,Y} = ^r ] (f). 

Thus, 

Y(K) = ±Y(7j(f)) = ±\y,fj](f)=0 
since [Y,77] ft = because 7i*[Y,fj] h = [n^Y, 77] = [0,77]. Hence, 

n // cy i = ^ a ij {V Di Dj, fj) + n = J2 <7 ij '<V Vi Vj, 77) + k = (n - l)if r + « 

where ifr is the mean curvature of T in M. 

In the sequel, we deduce some useful properties of the distance function 
d = dist(-,K) from K. We denote by Y e and K e the level sets d — e in 
M and M, respectively. Thus, T e and fT e are equidistant from Y and K, 
respectively. It is immediate that K e is a Killing cylinder over r e . Since T is 
assumed to be C 2 ' a , the function d is also C 2,a at points of \I/(R x fi e ), where 
f2 e C f2 is a tubular e-neighborhood of T in M for small e > 0. 

Given q e \1/(R x Q e ) we write q = exp p drj for some p & K. Hence, 



|£>i| = |vd| = 1. 



It follows that 



= ^A(W, W) = (V A Vd, ^'^-(d)^) = d k di, k . 
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We also have 



(V Do D ,Vd) := - Ke 



and 




Therefore, 



Ad 



-K e - a 1 " b\i = -nHi 



(13) 



where by are the components of the Weingarten operator A e and the 
mean curvature of K e . 

Fact 2 . All of the above calculations on the distance function d remain valid 
if we replace fl t by the larger subset Qq in Q consisting of the points which 
can be joined to T by a unique minimizing geodesic. It was shown in [12] 
that in this set d has the same regularity as T. 

In this paper the ambient Ricci tensor in direction v is defined by 



where R is the curvature tensor in M and ei, . . . ,e n ,v is an orthonormal 
basis. We follow [H] or [15] and use the result in Fact [2] for the proof of the 
following result. 

Lemma 3. Assume that the Ricci curvature satisfies Ric^- > — ninfr H^. 
Let y G T be the closest point to a given point x E T t C fi . If H cy \ > 0, 
then, we have 



Proof: At d = e and since A is the unit speed of a geodesic, on one hand 
we have that 



n 




1=1 




(A t Di, Dj) = A(V A A, A) = (V Dl V fl ,Di,^) + (V Dj Di, V Dl Dj) 

= -i? e (A, A) - (V [0i , Dl] A, A) + (V Ci A, V Dl A) 
= -R t {D u DA - (V c . A, [A, A]) + (V A A, V Dl Dj) 
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where R e = (R(-, Di)D 1: -)\d=e- On the other hand, 
— (A e Di,Dj) = {V D ^A e D i ,D j ) + {A t D i ,VD 1 D j ) 

= (A' e D h Dj) - {V Di D u V Dl Dj) - {V Dj D x , V Dl A). 
Adding the above equations we have the Ricatti equation 

A' e ~A 2 -R e = 0, 

Taking traces, we obtain 

n^-H^ = Di(tr A e ) = tiW Dl A e = tv(A 2 e + R e ) > n(H* cyl ) 2 + Ric#pi). 

From our hypothesis on Ric^- we have that z(d) = H cy \(d) — if cyl (?/o) satisfies 
z'(d) > H 2 yl (d) - mfH 2 cyl > H 2 yl (d) - H 2 cyl (y ) = (H cyl (d) + H cyl (y ))z(d). 

Since if cy i > 0, it follows that z'(d) > cz(d) in some interval d G [0, do > 0] 
for a constant c > 0. We obtain that Hl yl does not decrease with increasing d. 
This concludes the proof of the lemma. 



4 The C u estimate 

In this section, we obtain apriori C° estimates for solutions of the Dirichlet 
problem 

We construct barriers for u in ([8]) on Q (see Fact [2]) by 

ip(x) = sup + h(d(x)) 
r 

where d = dist( • , T) is regarded as the distance from T on M and the function 
h will be chosen later. We work with the frame V! := Vd, v 2 , . . . , v n and the 
corresponding frame D , D\, . . . , D n . Thus, 

Di(d) = (Vd, Di) = (D 1 ,D i ) = (vi,Vi) = (Vd,Vi) = Vi{d). 

We have, 

(Pi = h'di and = h"didj + h'd^j. 
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We obtain from (gj) and ([II]) that 

$ = - Si ) = ^(Zi'dj - Si) = h'<T - s i 

and 

<^iy = Viy - s *;j + r 7ij = hf'dkdj + h'di.j - Si-j + - 7y. 
Since 7^ is skew-symmetric, we have from (j3J) and that 

= W 2 (/i" + ti&i - 4) - - s*) - s?){h"didj + - s i;i ) 
= (W 2 - h' 2 + 2h'(Vd, Vs) - (Vd, Vs) 2 ) fc" + W 2 h'd] l + i2 

where 

R ■= -tij^drj - W 2 s i . i + (h a d*d j - h'{d i s j + dV) + sV><y. (14) 
Using that 

V^ 2 = / + £*£ fc = fj r h n - 2h'(Vd, Vs) + |Vs| 2 , (15) 
we conclude that 

= (/ + I Vs| 2 - (Vd, Vs) 2 ) /i" + W 2 h% + R (16) 

where R is a polynomial of second degree in h' and its coefficients are just 
functions on M. 

We have from (TT5]) that 

4 = <r ij d i;j = a"(V/, Vr/. Dj) = ^ = Ke - n#^. 

Since Dip = 7r*(/i' Vd — Vs) from (j4l , we also have 

(7r*V Do A), £>¥>) = /i'« e - (Va, A), Vs). (17) 

Thus, we obtain 

W 3 Q[ip] = (/ + I Vs| 2 - (Vd, Vs) 2 )//' - (/« e + n^ 2 # c e yl )/i' + if 

where 

ir :=i? + (/ + w 2 )(v Do A),vs). 
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We choose for PJ the test function 



p ca 

where A > diam(f2) and C > is a constant to be chosen later. Then, 

h ' = e C(A-d) and h »=_ Ch >. 

Hence, 

Observe that f/W 2 < 1. Moreover, as C — > oo we have that 1/W — > and 



1. 



W (/i' 2 -2/i'(Vrf,Vs) + |Vs| 2 + /) 1 /2 
In particular, we have 

— — ^ as C — > oo. 

Choose C ^> such that, in particular, C+/t e > 0. Using sup^ \H\ < infr i? cy i 
and Lemma [3l we obtain 

Q[y?] < -n\H\ < nH. 

Thus, one has at points of Qq that 

Q[f] < Q[u] = nH, (f\ r > u\r- 

We now prove that (p > u on Q. By contradiction, assume that there exist 
points for which the continuous function u* := u — (p satisfies u* > 0. Hence 
m := u*(y) > at a maximum point y G H of w*. Choose a minimizing 
geodesic 7 joining y to T for which the distance d = d(y, T) is attained. 
Thus, j(t) = exp yo trj, < t < d, starts from a point yo G T with unit 
speed T]. Since 7 is minimizing, we have d(j(t),T) = t and the function ip 
restricted to 7 is differentiable with (p'(j(t)) = e c ( A ~ l \ Since the maximum 
of u* restricted to 7 occurs at t — d, i.e., at the point y, one has that 

u'(j(d)) - cp'(j(d)) = (uj(j(d)) > 0. 
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This implies that 

(Vu(y), 1 '(d))>p'( 1 (d)) = e c ^>0. 

In particular Vu(y) ^ 0, and hence the level hypersurface 

S = {x G Q n B r (y) : u(x) = u(y)} 

is regular for small radius r. Along 5 we have 

u*(x) + <p(x) = u*{y) + <p(y) > u*(x) + <p(y), 

and since ip is an increasing function of d( ■ , T) we have d(x, V) > d(y, V) = d. 
From this we conclude that the points in S are at a distance at least d from 
T. Since S is C 2 it satisfies the interior sphere condition: there exists a small 
ball B e (z) touching S at y contained in the side to which Vu{y) and 7 ; (eQ 
points. Thus, the points of B E (z) satisfy u(x) > u(y), and hence 

<p(x) + m > u{x) > u{y) = (p(y) +m, x G B e (z), 

where in the first inequality we used the definition of m. Again because p 
is an increasing function of d, we have d(x,T) > d on B e (z) and therefore 
this ball is contained in the interior of Q far away from T. This allows us to 
extend the geodesic 7 through B £ (z). We claim that the center z of the ball is 
contained in this extension. Otherwise, the broken line consisting of 7 and of 
the radius in B £ (z) from z toy has length smaller than a minimizing geodesic 
joining z to y G T (for a suitable small e such a geodesic must cross the level 
hypersurface S at a point x 7^ y at distance to T greater than d). Thus, if 
there exists at least two distinct minimizing geodesies joining y to T, then the 
point z is contained in the extension of both geodesies after its intersection 
at y. Choosing e sufficiently small, we see that this configuration is not 
possible (the construction we made above applies to both geodesies). This 
contradiction implies that the maximum point y belongs to f2 . However, in 
this case, u*(y) < 0, a contradiction. We conclude that u < <p throughout fl 
and therefore ip is a continuous super-solution for the Dirichlet problem (jSJ). 

In a similar way, we may construct lower barriers for u, that is, continuous 
sub-solutions for (jSJ). It is a standard fact that the existence of these barriers 
implies the desired C° apriori estimates. Thus, we have proved the following 
result. 
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Lemma 4. Under the assumptions of Theorem^ there exists a constant 
C = C(n, H) such that 

Mo < C+ \<f>\o 
ifue C 2 (Q) n C°(0) satisfies Q[u] = nH and u\ T = 0. 



5 Boundary gradient estimates 

In this section our task is to produce apriori gradient estimates along T for 
the Dirichlet problem ([8]). This is accomplished by constructing local lower 
and upper barriers for E in a tubular neighborhood of T. 

We construct barriers of the form w + along a tubular neighborhood Q e 
of r as defined in Section 3. Here, w = ip(d(x)) for some real function ip to 
be chosen and d = dist( • , T). Moreover, the boundary data is extended to 
a function in Q e along the normal geodesies in a way we make precise later. 

We denote 

Q[u] = Q[u] - nH. 
A simple estimate using (fT2|) and then (|T7|) gives 

Q[w + 0] = a ij {x, Vw + V(f))(w i;j + 4> i;j ) + b(x, Vw + V0) - nH 

< a ij w iv + ^\^\ 2 , a + b-nH, (18) 



where 



and 



^ ■■=§-, = w aij - + + ^ (19) 

f + w 2 

(^'« e + (7T*V Do D , D<j>) - (Vd„ A), Vs)) 



W 3 

since n t = (Vd Dq, Vd) and 

D(w + 0) = cr^(^'^ + 0j - Sj)vj = D<p + 7r,(^' W - Vs). 

From now on Rj, j > 1, denotes a polynomial of at most second degree in 
ip' whose coefficients are functions in M. As in (ITB1) and ( 1T61) we first obtain, 

iy 2 = / + tjj' 2 - 2t/j'(Vd, Vs - V0) + | Vs - V0| 2 , 
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and then 

W 2 w l . t - w l w J w VJ = (f + \Vs - V0| 2 - (W, Vs - V0) 2 )^" + + Ri- 

Moreover, 

w l (f) j Wi-j = {i>'d l - s l )(i) j (4>" didj + ip'di-j - Si-j) 

= -V»V"<Vd, Vs - V0) + V"<VtZ, Vs) (W, Vs - V0) + R 2 

and 

ft&Wi.j = (f) l (f) j (ilj"did j + ip'di-j - Si-j) + R 3 . 

Now define 

i[>(d) =fi\n(l + Kd) 
for constants /x > and X > to be chosen later. We have 

V'=t^ and V" = -V- 
1 + a a /x 

Then using d^ = —nH^ yl + K e we obtain 

W 2 w l . t - ic'ir'ic,., = -ij'(nH e cjl - n t )W 2 + R 4 , 



and 

ft&Wi-j = R 6 . 

Since ( TT9T) gives 

W^a %i Wi-j = W 2 w % ;i - w l w j Wi-j - (w l (f) j + w j (p^ibi-j - ftfowij, 
we now conclude from (1TB]) that 

^ 3 Q[w + 0] < -^'(n^ c £ vl -K £ )^ 2 --^' 3 (W, Vs-V0) + (6-rxi/)W 3 + J R 7 . 

/x 

From the expressions above for b and W 2 it follows that 

tW 3 + tp'K € W 2 = R s . 
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Hence, we obtain 

2 



W 6 Q[w + 0] < -(n(H + H e cyl ) + -(W, Vs - V0))^ + R 9 . 

We choose \l in such a way that /i — > as K — > oo. Namely, 

C 



ln(l + AT) 

for some constant C > to be chosen later. As K —>■ oo we have that 

It also holds that ip'/W ~ 1 as K — ► oo. Thus, at points of T the last 
inequality becomes 

lV 3 Q[w + 0] < -(n(H + H cyl ) + -(Vs-V0,ri))ij ,3 + R 9 . 
We choose the extension of <fi in such a way that at points of V it holds 

(V0,7?> < (VS,77). 

Therefore, assuming that if cy i + if > and choosing K large enough, we 
assure that Q[w + <p] < on a small tubular neighborhood Q e of T. Notice 
that (w + 4>)\r = <f>\r- Choosing C and K large enough we also have that 
w + > w|r e + 0- Therefore, w + is a locally defined upper barrier for the 
Dirichlet problem (jSJ). A lower barrier may be constructed in a similar way. 
Thus, we have proved the following fact. 

Lemma 5. Assume that u e C 2 (Q) H C 1 (f2) satisfies Q[u] = nH and 
u\y = 0. If \u\ is bounded in fl, then 

sup | Viz | < C 
r 

by a constant that depends on \u\q. 
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6 Interior gradient estimates 



6.1 The prescribed mean curvature case 

In this general case, we adopt ideas from the classical estimate of Korevaar 
[T3] . Suppose that the maximum of \Du\ is attained at an interior point, 
say xq G fi, where we may assume that \Du\ 7^ without loss of generality. 
Consider a geodesic ball B = B(x ,p) C Q centered at x with small radius 
p < 1 so that I Dm I > C at points of B for some positive constant C . Without 
loss of generality, we may assume after a translation along the flow lines of 
Y, if necessary, that u < at points of the solid cylinder 7r _1 (.B). 

Let rj(x, s) > be a continuous function defined infix IR~ that vanishes 
on dB x IR~ and is smooth wherever it is positive. Then, let E be the normal 
geodesic graph over E defined by 

q = exp p erj(p)N(p) 

where p G E is parametrized by (x,u(x)). Recall that N given in (jSJ) was 
fixed to be upwards. 

For small e > 0, we may describe E as a Killing graph of some function u 
defined in Q. We denote by y the point in Q that maximizes u — u. It is clear 
that y G B and that DiU = DiU at this point. From ([3]) the tangent planes 
to both graphs have the same slope with respect to the flow line 7r _1 (y) of Y . 

We claim that 

Hu{y) < H u (y) (20) 

where H u and H u denote the mean curvature of E and E, respectively. In 
fact, moving E upward along the flow lines until the points (y,u(y)) G E 
and (y,u(y)) G E coincide, we obtain a tangency point for both graphs. 
Moreover, by the choice of y it is clear that the translated copy of E is above 
E locally around the point. Thus, the inequality (120]) is consequence of the 
comparison principle for the mean curvature PDE. 

In analytical terms, the above geometrical reasoning is justified in the 
following way: one has H u (x) = Q[u](x) and H u (x) = Q[v](x) since both 
hypersurfaces are described as Killing graphs. By construction, u = u at OB 
and u < u in B (this time, we are not considering the translation of the geo- 
metric proof). If Q[u] > Q[u] in B, then the analytical comparison principle 
(cf. Thm. 10.1 em [11]) assures that u < u in B. Thus, this contradiction 
shows that (|2"U|) holds. 
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It is a well-known fact that since the variation of £ we consider is along 
the normal direction, then the mean curvature may be expanded as 

nH u (x) = nH u (x) + eJr/ + 0(e 2 ), (21) 

where (x,u(x)) and (x,u(x)) parametrize correspondent points in E and £ 
along the same normal geodesic and 

J = A^ + \A\ 2 + RicM(N,N) 

is the Jacobi operator produced by the linearization of the mean curvature 
equation. Here, A s is the Laplace-Beltrami operator on £ and \A\ denotes 
the norm of its second fundamental form. 

Let x = y for some x. It follows from (1201) and (1211) that 

eJ V + 0(e 2 ) = n(H a (y) - H u (x)) < n(H u (y) - H u (x)). 

On the other hand, Taylor's expansion of H u gives 

H u (y) = H u (x) + e V H i T i + 0(e 2 ), 

where T % are the components of the horizontal projection of the normal vector 
field N. Thus, we get at y that 

A E r/ + \A\ 2 rj + Ric(iV, N)r} < nrjHiT + 0(e). 

Therefore, 

A s r? - M V < 0(e) (22) 

for some constant M > which does not depend on rj. 

Next we proceed as in [T3] choosing rj = g(6(x, s)) for some real function 
g to be chosen and a function 9 defined so that A^?7 is large for sufficiently 
large \Du(x)\. Since e is chosen small, then (|22|) will give a contradiction. 
Observe that C being large implies that the tangent hyperplanes to £ near 
(y,u(y)) are very sloppy. 

That a tangent hyperplane to £ is almost vertical means the tangential 
component of the gradient of 9 is approximately 9 S . Then, we define 

9(x,s) = (Ks + (p 2 -r 2 ))+ 

for some small constant K > 0, where r(x) = dist^^o, ^) is the geodesic 
distance measured from the center xq of B and ( • ) + m eans positive part. We 
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have that < 9 < p. Since we are assuming height estimates for £, we may 
choose K sufficiently small in such a way that 9 > in a neighborhood of 
(y,u(y)) in B x Mr . We restrict ourselves to points where 9 is differentiable. 
There, 

9 S = K > 0. 

Since 

A E 77 = </'|V s 0| 2 + </A s 0, (23) 
we have from (1221 and (1231) that 

#"|V s fl| 2 + g'Ax9 -Mg< 0(e). (24) 

By hypothesis, the tangent plane of £ at (y,u(y)) is not horizontal. Other- 
wise, we obtain from ([3]) that Du{y) = 0. Let e be the unit vector that gives 
the steepest ascent direction in the tangent plane of £ at (y,u(y)), namely, 

i\Du\ 2 D + fW&Dj). 



W\Du\ 

Denoting by V6 1 the ambient gradient of 9 and using that p < 1, we have 

f 1/2 / v j D(9)\ f 1/2 

<V E 0, e) = (W, e) = J w (K\Du\ + > ^(KW -CK-2), 

where C > is a constant independent of ii given by the following estimate: 



II: 



Dj(9) = ——(KDjis) - 2r Vj (r)) > -2 - CK. 



\Du\ Jy ' \Du 

Since K and C are independent of u and the parameter s, we may assume 
that | Du | > 2/K + C, and conclude that 

|V E 0| > 0. 

Finally, for Ci > large we choose 

g(9) = - 1. 

It is easily seen that this choice leads to a contradiction with (|24p . We 
conclude that |£)tt| and therefore |Vm| is bounded by some constant which 
does not depend on u. 
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Lemma 6. Assume that u G C 3 (Q) D C 1 (f2) satisfies Q[u] = nH and 
u\y = 0. Ifuis bounded in Q and |Vit| bounded in T, then \ Vu\ is bounded 
in Q by a constant that depends only on \u\o and sup r |Vit|. 

The usual elliptic regularity results guarantee that the above estimate is 
also true for a C 2,a function (see 1111). 



6.2 The constant mean curvature case 

In this case a standard argument works. In fact, consider the positive function 

/ 



6:= (N,Y) 



y r FT\D 



u\ 



since a lower estimate of G clearly yields an upper estimate for \Vu\. Under 
the assumption that H is constant and being Y a Killing field, it is well- 
known (cf. |2j) that G is a Jacobi field, namely, JG = 0. By assumption 
the Ricci tensor is bounded from below. Thus, since S is compact there is a 
constant c > such that \A\ 2 + Ric A 7/(iV, N) > —c. Thus G is a supersolution 
to the elliptic operator As — c. Hence, the classical minimum principle states 
that 

min Q > min G. 

s as 

This assures that |Vw| is uniformly bounded from above by a constant in- 
volving the boundary estimates for |Vtt|. 



7 The proof of the theorem 

In view of the Continuity Method, one must seek for an initial minimal 
surface with boundary given by T. This may be accomplished by defining 
the sets 

C = {u £ C 0,1 (il) :u\ r = 0} 

and, given k > 0, 

Cfc = {u G C : \u\o,i < k}. 

The hypothesis on the existence of an immersion l: Q — ► M assures that 
the set C is non-empty since we may consider the hypersurface as the 
graph So of the function u — 0. For the = 0, if we assume M 
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is geodesically complete, the immersion i may be obtained as follows: we 
construct a geodesic cone by joining points of a Killing graph in K over V 
to a vertex po inside Mo. This cone is contained inside Md since the Killing 
cylinder K is mean convex and Md is geodesically complete. Moreover, it 
may be smoothed out near the vertex. The resulting hypersurface may be 
given as a Killing graph since the geodesic cone is always transversal to the 
geodesic vertical fibers. 

We then formulate the issue of the existence of a minimal graph spanning 
i(T) as the minimization of the functional 

T(u) = / W{x, Vu(x)) v/adx, ueC, 
Jn 

where 

W = v 7 / + UiU* 

and the first derivatives of u are taken in a weak sense. Notice that / and 
Ui = Ui — f x l 2 bi do not depend on u. It is clear that u is a critical point of X 
if and only if is a weak solution of the mean curvature equation in divergence 
form. Since the principal part of the mean curvature equation is positive- 
definite and the coefficients of this equation (including H) do not depend on 
the function, it follows from Theorem 11.10 and Theorem 11.11 in [11] that 
X has a extremum in C. In fact, these theorems require upper bounds in 
the Lipschitz norm of the candidates mgC which may be obtained from the 
apriori C 1 estimates we derived earlier. 

The C 2 ' a regularity of the minimizer function uq follows from very general 
results found in [T3]. This function defines a minimal graph over Q with 
boundary T. 

For the proof of the existence part we apply the well-known continuity 
method to the family of Dirichlet problems 

Q a [u] = naH, u\ r = cr0, 

where a G [0, 1]. The subset / of [0, 1] consisting of values of a for which there 
is a solution is non-empty since we have an initial minimal graph spanning 
the boundary data (f>. The openness of J is a direct consequence of a standard 
application of the implicit function theorem since the derivative of Q a is a 
linear homeomorphism. The closedness of I follows from the apriori estimates 
we had proved and linear elliptic PDE theory. Thus, the continuity method 
assures that lei. 
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In order to prove the uniqueness statement, we deduce a kind of flux 
formula. We suppose that there exists a hypersurface E' in M with <9E' = T 
and whose mean curvature is the same as E at corresponding points in flow 
lines. This means that if x = ir(p) for p e S' then the mean curvature of E' 
at p is H(x). Translating E' we may suppose that E, E' and a part of the 
cylinder K form an oriented cycle which bounds a domain U in M . Since 
Y is tangent to the part of the boundary of U contained in K, we conclude 
from divergence theorem applied to the field BY in U that 



where N and N' define respectively the orientations in E and E'. Applying 
now the divergence theorem to the hypersurfaces E and E' we obtain that 



where v and v' are respectively the outward unit co-normals to T with respect 
to E and E'. This implies that there exists a point p in T where E and E' 
are tangent, that is, where v\ p = v'\ p . Thus, since E' is locally a graph near 
p, we conclude from the the boundary maximum principle that E = E'. This 
concludes the proof of the theorem. 
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